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Local  skin  friction  reductions  have  been  measured  using;  an  array  of 
flush-mounted  hot-film  probes  in  a  microbubble-mod  if ied ,  zero  pressure 
gradient,  turbulent  boundary  layer.  The  results  of  earlier  integrated 
skin  friction  measurements,  that  showed  the  reduction  to  be  a  function 
of  plate  orientation,  gas  flow  rate  and  f  roes t ream  velocity,  have  been 
confirmed  both  qualitatively  and  quantitatively.  With  the  measurement 
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for  as  much  as  60-70  boundary  layer  thicknesses  downstream  of  the  injection 
region.  It  is  further  shown,  using  a  probe  flush  mounted  just  upstream  of 
the  injection  section,  that  there  is  no  apparent  upstream  interference  due 
to  the  gas  injection.  Spectral  measurements  indicate  that  microbubbles 
can  cause  a  reduction  of  high  frequency  shear-stress  fluctuations.  This 
suggests  a  destruction  of  some  of  the  turbulence  in  the  near  wall  region. 
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From:  N.  K.  Madavan,  S.  Deutsch  and  C.  L.  Merkle 

Subject:  Measurements  of  Local  Skin  Friction  in  a  Microbubble  Modified 

Turbulent  Boundary  Layer 

Abstract:  Local  skin  friction  reductions  have  been  measured  using  an 

array  of  flush-mounted  hot-film  probes  in  a  microbubble-modified,  zero 
pressure  gradient,  turbulent  boundary  layer.  The  results  of  earlier 
integrated  skin  friction  measurements,  that  showed  the  reduction  to  be  a 
function  of  plate  orientation,  gas  flow  rate  and  freestream  velocity,  have 
been  confirmed  both  qualitatively  and  quantitatively.  With  the  measurement 
plate  above  the  boundary  layer,  it  is  shown  that  skin  friction  is  reduced 
monotonical  Iv  for  all  air  flow  rates  at  each  of  three  freestream  velocities 
between  4  and  17  m/ sec .  For  the  plate  below  the  boundary  layer,  however,  it 
is  possible  for  increasing  gas  injection  to  lead  to  smaller  local  skin 
friction  reduction  at  the  lowest  speeds.  Drag  reduction  appears  to  persist 
tor  as  much  as  60-70  boundary  layer  thicknesses  downstream  of  the  injection 
region.  It  is  further  shown,  using  a  probe  flush  mounted  just  upstream  of  the 
injection  section,  that  there  is  no  apparent  upstream  interference  due  to  the 
gas  injection.  Spectral  measurements  indicate  that  microbubbles  can  cause  a 
reduction  of  high  frequency  shear-stress  fluctuations.  This  suggests  a 
destruction  of  some  of  the  tiubuience  in  the  near  wall  region. 
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I.  INTRODUCTION 

The  authors  (Madavan,  Deutsch  and  Merkle  1984a,  1984b;  hereafter 
referred  to  aS  MDM1  and  MDM2)  have  recently  shown  that  the  introduction 
of  gas  bubbles  (microbubbles)  into  a  hydrodynamic  turbulent  boundary  layer 
results  in  significant  reductions  in  skin  friction.  In  a  series  of 
experiments  building  on  some  earlier  Soviet  work  (Migirenko  and  Evseev 
1974;  Bogdevich  and  Malyuga  1976;  Bogdevich  and  Evseev  1976),  integrated 
skin  friction  measurements  were  made  on  a  102  mm  x  254  mm  force  balance 
downstream  of  a  porous  injection  section.  Microbubbles  were  introduced 
into  the  boundary  layer  by  injecting  gas  (usually  air)  through  the  porous 
surface.  The  effect  of  buoyancy,  gas  flow  rate,  freestream  velocity, 
gravitational  orientation  and  porous  material  on  the  integrated  skin 
friction  reduction  have  been  documented  in  MDMl  and  MDM2.  Generally 
speaking,  the  measurements  have  shown  that  the  magnitude  of  the  skin 

friction  reduction  increases  with  air  flow  rate  and  that  it  is  the  ratio 

of  the  volumetric  flow  rate  of  gas  to  that  of  water  that  is  the 

important  parameter.  Tests  with  various  types  of  porous  surfaces, 
including  pore  size  variations  over  a  range  of  0.5  to  100  pm,  have  >-'nown 
that  the  skin  friction  reduction  can  be  realized  with  a  variety  > '  porous 
materials,  although  details  of  the  variation  of  the  skin  friction 
reduction  with  air  flow  rate  are  slightly  dependent  on  the  type  of 
material . 

The  microbubble- laden  turbulent  boundary  layer  is  an  extremely 
complicated  flowfield  and  the  mechanisms  which  read  to  the  observed 
skin  friction  reductions  remain  unclear.  Certainly  a  proper  under¬ 
standing  of  the  phenomena  cannot  be  expected  from  integrated  skin 
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friction  measurements  alone  and  calls  for  the  measurement  of  additional 
details  of  the  flowfield.  A  more  precise  knowledge  of  the  raicrobubble 
volumetric  concentration  profiles  and  bubble  size  distributions,  as  well 
as  the  modified  mean  velocity  profiles  and  corresponding  turbulence 
characteristics,  are  necessary  prerequisites  for  achieving  such  under¬ 
standing.  This  information,  however,  is  very  difficult  to  obtain.  First, 
it  is  noted  that  the  bubble  concentrations  of  interest  are  quite  high 
(Migirenko  and  Evseev  1974  have  reported  peak  concentrations  as  high  as 
50  to  80%  by  volume)  so  that  the  microbubble-laden  boundary  layer  is 
optically  opaque,  thus  precluding  the  use  of  contemporary  optical 
measurement  techniques.  Second,  both  the  earlier  reported  gravitational 
effects  (MDM1 ,  MDM2)  and  simple  models  of  the  bubble  formation  process 
suggest  the  phenomenon  is  operative  only  at  high  unit  Reynolds  numbers. 
Thus,  the  boundary  layers  of  interest  are  extremely  thin,  further 
aggravating  measurement  attempts  (boundary  layer  thicknesses  are  less 
than  10  mm,  while  sublayer  thicknesses  are  a  few  tens  of  micrometers  in 
the  current  experiments). 

Measurements  of  bubble  concentration  profiles  reported  by  Migirenko 
and  Evseev  (1974),  as  well  as  the  authors'  own  results  (MDM1),  suggest  a 
bubble-free  layer  very  close  to  the  wall.  Local  skin  friction  measurements 
using  flush-mounted  hot-film  probes  become  then  a  practical  and  attractive 
first  step  toward  understanding  the  local  characteristics  of  the  micro¬ 
bubble  boundary  layer.  In  the  current  study,  an  array  of  flush-mounted 
hot-film  probes  has  been  used  to  investigate  the  downstream  evolution  and 
persistance  of  the  skin  friction  reduction.  The  instantaneous  signals 
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from  the  hot  films  have  also  been  utilized  to  compare  the  spectral  content  of 
the  turbulence  in  the  presence  and  absence  of  microbubbles. 


II.  EXPERIMENTAL  DETAILS 

(a)  Facility  and  Setup 

The  experiments  were  conducted  in  the  rectangular  test  section  of  the 
water  tunnel  facility  at  The  Pennsylvania  State  University  using  the  boundary 
layer  on  the  tunnel  test  section  wall.  The  tunnel  is  a  closed  circuit  system 
and  the  rectangular  test  section  has  dimensions  of  508  mm  x  714  mm  x  762  mm 
long.  The  flat  plate  on  which  the  measurements  were  performed  was  designed 
to  replace  an  existing  279  mm  x  533  mm  long  test  section  window.  A 
102  mm  x  178  mm  long  porous  section  for  gas  injection,  followed  by  a 
102  mm  x  254  mm  force  balance  was  centered  spanwise  on  this  test  plate 
(see  Fig.  1).  Details  of  the  facility  used  and  the  experimental  setup  can 
be  found  in  MDM1 . 

An  array  of  six  hot-film  probes  (TS1  147 1 W )  were  employed;  five  of 
these  were  flush  mounted  on  the  force  balance  downstream  of  the  injection 
section,  while  the  sixth  probe  was  mounted  immediately  upstream.  The 
locations  of  the  probes  are  shown  schematically  on  Fig.  2.  The  probes 
were  connected  to  constant  temperature  anemometers  (both  DISA  55M01  and 
TSI  1750  sets  were  used).  Normally,  the  probes  were  operated  at  an 
overheat  ratio  of  1.06,  this  value  being  chosen  as  a  compromise  between 
probe  sensitivity  and  probe  life. 

The  porous  material  used  for  air  injection  was  of  sintered  stainless 
steel,  commercially  available  for  use  as  filter  material.  Experiments  were 
conducted  with  the  smallest  available  filter  material  (0.5  um) ,  as  well  as 
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with  a  larger  size  (100  pm).  The  numbers  in  parentheses  refer  to  the  smallest 
particle  size  that  can  pass  through  the  filter,  not  the  actual  pore  size. 

The  tunnel  test  section  velocity  was  monitored  by  a  transducer  that 
measured  the  pressure  drop  across  the  tunnel  contraction  section.  Injected 
air  flow  rate  and  pressure  were  measured  using  a  flowmeter  and  a  pressure 
transducer . 

( b)  Data  Acquisition  and  Run  Procedure 

Two  methods  of  data  acquisition  were  used.  In  the  initial  experiments, 
the  signals  from  the  six  hot  films  were  recorded  on  three  Nicolet  digital 
oscilloscopes  and  stored  on  floppy  disks.  The  storage  capacity  of  the  disks 
limited  these  measurements  to  2048  points.  These  digitized  signals  were  then 
transferred  to  an  IBM  4341  mainframe  system  for  reduction  and  analysis.  Other 
flow  parameters,  namely,  tunnel  freestream  velocity,  injected  air  flow  rate 
and  pressure,  were  taken  manually  from  digital  voltmeter  readouts  at  the 
instant  the  hot-film  signals  were  recorded.  In  the  latter  stages  of  the 
experiments,  an  analog-to-digital  (A/D)  converter  was  used  to  digitize  the 
hot-film  signals  and  feed  the  data  directly  to  an  on-line  VAX  11/780  computer. 
Here,  the  allowable  sample  size  was  much  larger.  For  both  methods  of  data 
acquisition  the  sampling  rate  was  typically  chosen  as  2,000  Hz.  The  minimum  ■ 
record  length  used  was  about  1  second.  In  both  procedures,  the  data  channels 
were  sampled  sequentially,  the  maximum  delay  between  the  first  and  last 
channel  sampled  being  about  10  us.  The  in'stantaneous  voltages  from  the 
films  were  recorded  and  linearization  was  accomplished  on  the  computer. 
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The  run  procedure  used  for  the  local  skin  friction  measurements  reported 
here  was  as  follows.  The  tunnel  velocity  was  held  constant  and  the  air  flow 
rate  was  set  at  levels  ranging  from  zero  to  maximum.  At  each  air  flow 
setting,  the  oscilloscopes  were  simultaneously  triggered  (or  sampling  was 
initiated  on  the  A/D  converter)  to  record  the  hot-film  signals.  Other  flow 
parameters  were  recorded  at  the  same  time  the  oscilloscopes  were  triggered. 
(With  the  A/D  converter,  these  parameters  were  also  automatically  sampled.) 
Immediately  before  and  after  each  run  with  microbubbles,  probe  calibration 
runs  were  performed  by  recording  the  hot-film  signals  over  the  range  of 
tunnel  velocities  in  the  absence  of  microbubbles. 

( c)  The  Experimental  Boundary  Layer 

The  boundary  layer  into  which  microbubbles  were  introduced  was  documented 
by  measuring  the  mean  velocity  and  turbulence  intensity  profiles  and  stream- 
wise  pressure  gradient  in  the  tunnel  test  section.  Details  can  be  found  in 
MDM1  ,  but  for  the  sake  of  completeness  and  the  important  role  these  data  play 
in  hot  film  probe  calibration,  the  major  results  are  highlighted  here. 

Mean  velocity  profiles  were  obtained  at  three  streamwise  locations  (shown 
in  Fig.  2)  and  at  two  different  tunnel  speeds  (4.7  and  10.6  m/s)  using  a  Laser 
Doppler  Anemometer  ( LDA) .  These  results  are  presented  in  wall  variables  in 
Fig.  3,  taken  from  MDM1 .  The  friction  velocity  used  in  the  figure  was 
evaluated  by  fitting  the  near-wall  data  to  Coles  (1968)  law  of  the  wall.  The 
boundary  layer  was  checked  for  "full  development"  by  examining  the  velocity 
profiles  in  outer  variables,  following  Purtell  et  al.  (1981).  This  is  shown 
in  Fig.  4,  also  taken  from  MDM1 .  Similarity  in  the  profiles  is  noted  at 
both  tunnel  speeds.  These  velocity  measurements,  along  with  measurements  of 
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the  pressure  gradient  along  the  tunnel  test  section  (see  MDM1 ) ,  establish  the 

experimental  boundary  layer  as  being  reasonably  close  to  a  classical,  fully 

developed,  zero  pressure  gradient  turbulent  boundary  layer. 

The  calibration  procedure  adopted  here  for  the  hot  films  (as  discussed 

below)  necessitated  establishing  a  virtual  origin  for  the  experimental 

boundary  layer.  This  was  accomplished  using  the  boundary  layer  displacement 

and  momentum  thicknesses  (6*  and  0),  evaluated  from  the  mean  velocity  profiles 

discussed  above.  Since,  for  a  flat  plate  turbulent  boundary  layer,  both  6 

and  0  are  proportional  to  the  distance  from  the  virtual  origin  to  the  0.8 

power  (4/5),  the  latter  can  be  linearly  extrapolated  from  a  plot  of  6*^^  and 
5/4 

9  versus  streamwise  distance  (referenced  to  a  known  location  along  the 
boundary  Layer).  By  this  procedure,  the  location  of  the  virtual  origin  was 
established  180  mm  upstream  of  the  porous  injection  section,  as  shown  in 
Fig.  5. 


(d)  Skin  Friction  Probe  Calibration  Procedure 

Flush-mounted  hot-film  probes  relate  the  heat  transfer  rate  from  the 

fiLm  to  the  fluid  shear  stress  at  the  wall.  For  constant  temperature 

operation  of  the  film,  this  relationship  between  the  instantaneous  voltage 

output  by  the  anemometer,  E  (which  is  directly  related  to  the  electrical 

heat  dissipated  bv  the  film),  and  the  instantaneous  shear  stress,  t  ,  is 

w 

of  the  form. 


.1/3 


AE2  + 


w 

Here  A  and  B  are  constants  that  must  be  determined  by  a  calibration 
in  .i  flow  where  the  shear  stress  is  known  either  from  theory  or  an 
measurement.  Brown  (19b7),  Liepmann  and  Skinner  (1954),  and  others 


(1) 

procedure 

independent 

have 


suggested  that  the  constants,  A  and  B,  determined  from  a  calibration  in  a 
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laminar  flow  can  be  used  for  measurement  in  a  turbulent  flow.  The  validity  of 
such  a  procedure  has,  however,  been  questioned  (see,  for  example,  Bellhouse 
and  Schultz  1966;  Bradshaw  and  Gregory  1961),  since  the  response  of  the  hot 
film  in  a  turbulent  flow  is  not  necessarily  identical  to  that  in  a  laminar 
flow  with  the  same  shear  stress.  Moreover,  it  is  difficult  to  obtain  a 
laminar  flow  for  the  calibration  experiment,  especially  one  that  can  cover  the 
range  of  shear  stress  values  expected  in  an  actual  turbulent  flow  application 
(Ramaprian  and  Tu  1983).  This  is  particularly  true  for  the  high  unit  Reynolds 
number  boundary  layers  of  interest  here. 

In  the  present  study,  an  in-situ  calibration  procedure  based  upon  the 
verified  characteristics  of  the  flat  plate  turbulent  boundary  layer  was  used. 
This  was  accomplished  by  recording  the  voltage  signals  from  each  probe  at  a 
series  of  different  tunnel  freestream  velocities  in  the  absence  of  micro¬ 
bubbles.  Each  signal  record  was  digitized  and  stored  and  the  time-averaged 

voltage  E  was  computed.  The  mean  local  skin  friction  value  r  at  each  probe 

w 

location  and  each  tunnel  velocity  was  then  obtained  from  classical  turbulent 
boundary  layer  correlations  of  skin  friction  and  length  Reyn'  ' ds  number,  using 
the  virtual  origin  established  from  the  LDA  measurements  of  the  velocity 
profiles.  The  use  of  the  classical  correlation  is  justified  by  the  velocity 
profile  measurements  that  indicate  the  present  boundary  layer  is  a  zero 
pressure  gradient  turbulent  boundary  layer.  The  specific  correlation  used 
as  the  one  suggested  by  White  (1974), 


w 


i/2rir'; 


0.455 


In  0.06  Re  ! 

x ' 


(2) 


where 


is 


the  densitv  of  water  and  Re  the 
'  x 


Reynolds  number  based  on 


the 


strenmwise  length,  x,  measured  from  the  virtual  origin  to  the  probe  location. 
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It  is  worth  noting  that  at  these  Reynolds  numbers  (roughly  2  to  10  million) , 
a  10%  measurement  error  in  virtual  origin  will  result  in  no  more  than  a  1% 
discrepancy  in  local  shear  stress. 

Having  obtained  this  tabular  set  of  mean  voltages  and  corresponding 
values  of  skin  friction,  the  calibration  constants  A  and  B  in  Eq .  (1)  are 
evaluated  by  an  iterative  technique.  In  this  procedure,  A  and  B  are  initially 
approximated  by  fitting  Eq .  (1)  to  the  time-averaged  tabular  data.  These 
values  are  then  iteratively  improved  by  the  method  described  by  Sandborn 
( 1 979a ; 1 979b)  and  Hanratty  and  Campbell  (1982).  In  the  present  study,  one 
iteration  was  found  sufficient.  Two  calibration  curves  were  fitted  to  the 
data,  one  at  low  and  the  other  at  high  velocities.  A  typical  calibration 
curve  is  shown  in  rig.  6. 

It  should  be  emphasized  here  that  the  calibration  measurements  were 

made  immediately  prior  to  and  after  measurements  with  air  injection  (sec 

Fig.  6).  The  average  of  these  two  calibrations  was  used  for  reducing  the 

hot-film  data  in  the  presence  of  microbubbles.  In  cases  where  the  pre- 

and  post-calibration  showed  deviations  which  resulted  in  the  correlation 

coefficient  of  the  linear  least-squares  fit  to  the  data  being  less  than 

0.95,  the  entire  data  set  was  discarded. 

For  the  case  when  there  are  no  bubbles  in  the  boundary  layer  a  measure 

of  the  vaLidity  of  the  hot-film  technique  can  be  obtained  by  comparing  the 

rms  intensity  of  the  skin  friction  fluctuations,  t* ,  measured  here  to  those 

w 

reported  in  the  literature.  Mitchell  and  Hanratty  (1966)  and  Sirkar  and 

Uanattv  (1970)  report  values  of  t'/t  ranging  between  0.32-0.36  as  obtained 

w  w 

from  electrochemical  shear  stress  measurements .  Kckelmann's  (  1974) 
measurements  in  the  thick  viscous  sublayer  of  a  glycerine  boundary  layer 
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using  a  flush-mounted  hot  film  provide  a  value  of  0.24.  By  assuming  that  very 
near  the  wall  0  <  y  <  0.1  }  both  the  streamwise  velocity  fluctuations  u'  and 
the  mean  velocity  U  vanish  linearly  with  decreasing  distance  from  the  wall, 
EckeLmann  noted  that  the  value  of  0.24  was  in  agreement  with  the  asymptotic 
slope  of  Lanier's  (1950)  measurements  of  u'/U. 

The  measured  value  of  t'/t  in  the  absence  of  microbubbles  are  compared 

w 

with  the  0.24  value  in  Fig.  7  for  the  tunnel  velocity  range  tested.  The  data 
were  typically  scattered  around  the  value  0.2.  The  agreement  with  earlier 
resuLts  seems  quite  reasonable,  particularly  when  one  considers  the  difficulty 
in  obtaining  resuLts  at  the  high  unit  Reynolds  numbers  and  thin  sublayers 
encountered  in  tins  studv. 

1  1 L .  KKSLLTS  AND  mSCLSSlON 

Experiments  were  conducted  at  tunnel  freestream  velocities  ranging 
between  4  and  17  m/ s .  For  the  most  part,  the  porous  section  fabricated  from, 
the  smallest  available  filter  material  (0.5  pm)  was  used.  Additional  experi¬ 
ments  were  also  conducted  using  the  100  pm  porous  material.  Air  flow  rate's 
ranged  up  to  0.(>06  m^/s.  As  is  shown  later,  this  flow  rate  corresponds  to 
a  volumetric  fraction  in  the  boundary  Layer  of  about  0.5,  assuming  that  all 
ot  the  air  ? 1 ov  remains  in  the  boundary  layer,  which  may  not  be  the  case. 

The  1  >eai  skin  friction  measurements  with  microbubbles  reported  here 

are  presented  in  terms  ot  the  skin  friction  coefficient,  C ,  normalized  to 

f 

the  correspond  n.g  skin  triction  coefficient  without  microbubbles ,  Cl  ,  as  \ 
f  uuvt ion  ot  the  volumetric  t  ruction  or  iir  in  the  boundary  layer.  The  air 
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flow  traction  was  defined  as  the  volumetric  flow  rate  of  air  divided  by  the 
total  flow  in  the  boundary  layer 


q 

a 

q  +  q 

a  w 

where 

q  =  u  f s  -  s* } b 

w 


(3) 


(4) 


and  the  subscripts  a  and  w  denote  air  and  water.  The  boundary  layer  thickness 
and  displacement  thickness,  6  and  5  ,  are  obtained  from  classical  correlations 
(see  White  1974).  These  thicknesses  were  calculated  at  the  mid-point  of  the 
porous  section  (in  the  absence  of  microbubbles).  The  parameter  b  is  the  width 
of  the  porous  section. 


( a )  Plate  Below  the  Boundary  Laver 

The  local  skin  friction  measurements  at  the  •>:  s  hot-film  probe 

locations  with  the  piate  below  the  boundary  layer  au  shown  in  Figs.  8 
through  10.  Note  that  this  orientation  (referred  to  as  plate  on  bottom  on 
the  figures),  corresponds  to  the  view  shown  in  Fig.  1,  where  the  experimental 
setup  is  mounted  on  the  bottom  wall  of  the  tunnel.  The  measurements  shown 
in  Fig.  8  are  for  a  tunnel  ve loci  tv,  U  =  16.7  m/ s ,  while  those  in  Figs.  9 

CO 

and  10  correspond  to  =  10.8  m/s  and  Ua  =  4.8  m/s,  respectively.  The 
figures  represent  measurements  using  the  0.5  pm  porous  section.  Skin 
friction  reductions  are  noted  at  all  tunnel  velocities,  but  the 
characteristic  trends  with  volume  fraction  depend  upon  IJ  .  The  reductions 

to 

are  greatest  in  the  region  close  to  the  porous  section  (probes  2  and  3); 
farther  downstream,  the  skin  triction  relaxes  hack  toward  its  undisturbed 
(no  mi c robuhb l es )  value  (probe  8).  This  is  probably  due  to  the  decreasing 
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air  content  in  the  boundary  layer  resulting  from  both  boundary  layer  growth 
and  bubble  migration.  At  the  high  and  intermediate  tunnel  velocities  (see 
Figs.  8  and  9),  the  decrease  in  skin  friction  is  seen  to  be  monotonic  with 
increasing  air  flow  except  possibly  for  probes  5  and  6  at  low  air  flows. 

In  contrast  to  these  high  and  intermediate  velocity  results,  the  low 

velocity  results  in  Fig.  10  show  that  the  skin  friction  goes  through  a  minimum 

and  then  increases  with  increasing  air  flow.  This  indicates  an  optimum  air 

flow  rate  for  maximum  skin  friction  reduction  at  low  velocities.  Based  on  the 

results  of  Silberman  (1957)  and  Hughes  et  al .  (1979)  that  indicate  the  bubble 

diameter  is  proportional  to  the  square  root  of  the  ratio  Q  /U  ,  one  can 

a  00 

speculate  that  the  reason  for  this  reduced  effectiveness  of  the  bubbles  at 
high  air  flow  rates  and  low  velocities  may  be  due  to  the  larger  bubble  sizes. 
These  larger  bubbles  quickly  migrate  out  of  the  boundary  layer  with  this  plate 
orientation,  reducing  their  effectiveness.  The  same  effect  was  not  observed 
at  the  high  velocities,  apparently  because  of  the  shorter  convection  time. 

( b)  Plate  Above  the  Boundary  Layer 

The  entire  test  section  of  the  water  tunnel  can  be  rotated  through  180 
degrees,  thus  effectively  positioning  the  experimental  setup  on  the  upper 
wall  of  the  tunnel.  At  this  orientation,  the  plate  is  above  the  boundary 
layer  (referred  to  as  plate  on  top  in  the  figures)  and  the  forces  due  to 
gravity  act  to  keep  the  bubbles  in  the  boundary  layer.  Earlier  reported 
force  balance  measurements  (MDM1)  showed  substantially  higher  reductions  in 
the  integrated  skin  friction  at  this  orientation,  especially  at  the  lower 
velocities.  The  hot  film  measurements  are  shown  in  Figs.  11  through  13  and 
compared  to  the  measurements  with  the  plate  below  the  boundary  layer  for  a 


-13- 


24  August  1984 
NKM: SD: CLM : Ihz 


range  of  tunnel  velocities.  Figure  11  compares  measurements  at  =  16.8  m/ s 
for  the  two  plate  orientations  at  probe  3  and  probe  6  (note  that  probe  3  i.-, 
located  close  to  the  injection  section  and  probe  6  is  located  farthest  down¬ 
stream).  While  the  skin  friction  reductions  at  probe  3  are  of  the  same 
magnitude,  slightly  larger  reductions  are  noted  at  the  downstream  probe 
location  when  the  plate  is  above  the  boundary  layer.  Figure  12  is  a  similar 
plot  at  =  10.7  m/ s  and  compares  the  measurements  for  the  two  orientations 
at  probe  2  and  probe  6.  Once  again,  it  is  noted  that  close  to  Lhe  injection 
section,  the  skin  friction  reductions  are  independent  of  gravitational 
orientation,  while  farther  downstream  the  plate-on-top  orientation  again 
results  in  somewhat  larger  reductions.  These  increased  reductions  with 
the  plate  on  top  occur  because  the  bubbles  remain  longer  in  the  boundary  layer 
where  they  are  most  effective.  Figure  13  presents  additional  data  for  a 
tunnel  velocity  at  4.b  m/s.  Note  that  data  for  probe  6  with  the  plate  above 
the  boundary  layer  are  not  available.  Quantitative  comparison  between  the  two 
orientations  in  Fig.  13  is  difficult  because  the  maximum  air  flow  rates  shown 
with  the  plate  above  the  boundary  layer  are  quite  small.  The  fluctuations  in 
the  hot-film  signals  at  higher  air  flow  rates  were  very  high  and,  due  to  the 
resulting  lack  of  confidence  in  the  measurments,  these  results  are  not 
presented.  The  data,  however,  suggested  that  the  skin  friction  reductions 
with  increasing  air  flow  were  monotonic,  as  compared  to  the  plate-on-bottom 
orientation  where  a  definite  "bucket"  in  the  curves  was  noted.  These  trends 
seem  reasonable  when  one  notes  that  the  maximum  integrated  skin  friction 
reductions  reported  in  MIJM1  were  observed  at  this  velocity,  =  4.6  m/s 

with  the  plate  above  the  boundary  layer. 
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( c)  Effect  of  Pore  Size  Variations 

The  size  of  the  microbubbles  is  certainly  an  important  parameter  in  the 
skin  friction  reduction.  Although  both  the  Soviet  results  (Bogdevich  and 
Evseev  1976)  and  intuition  suggest  that  the  pore  size  would  be  an  important 
variable,  literature  (Silberman  1957;  Hughes  et  al .  1979)  on  bubble  sizes 
in  turbulent  shear  flows,  though  sparse,  indicates  that  the  bubble  size  is 
not  determined  primarily  by  the  size  of  the  pores  used  for  injection,  but 
by  the  characteristics  of  the  flow.  To  investigate  this,  two  widely 
different  pore  sizes,  0.5  ym  and  100  ym,  were  used  for  gas  injection  to 
compare  the  magnitude  of  the  skin  friction  reductions.  Figures  14  and  15 
compare  measurements  made  using  the  100  ym  and  0.5  ym  porous  materials  for 
the  plate  below  the  boundary  Layer.  Figure  14  compares  data  at  probes  2 
and  6  at  UM  =  16.8  m/ s .  At  U =  10.7  m/s  a  direct  comparison  is  not 
available  and  for  clarity.  Fig.  15  presents  only  the  results  of  the  closest 
(probe  2)  and  the  farthest  (probe  6)  locations  for  the  0.5  ym  porous  section 
(presented  earlier  in  Fig.  9).  These  data  show  that  injection  pore  size 
does  not  have  a  substantial  effect  on  the  skin  friction  results  tor  the 
conditions  tested  here.  This  observation  is  in  agreement  with  the  integrated 
skin  friction  measurements  reported  earlier  in  MDM2.  These  results  are  of 
some  practical  importance  since  with  the  larger  pore  size  material  much  less 
energy  need  be  expended  to  introduce  a  given  quantitv  of  microbubbles. 

(d)  Persistence  of  the  Skin  Friction  Reduction  Phenomenon 

The  downstream  distances  to  which  the  skin  friction  reductions  persist 
was  studied  by  cross-plotting  the  data  reported  above  in  terms  of  a  non- 
dimensional  length  corrdinate,  x/5,  at  selected  air  flow  rates.  Here  x 
is  the  downstream  distance  measured  from  the  trailing  edge  of  the  porous 
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section  and  5  denotes  the  undisturbed  boundary  layer  thickness  at  this 

location.  These  results  are  shown  in  Figs.  16  through  19.  The  figures 

have  been  reconstructed  from  mean  lines  drawn  through  C'/C'  versus  Q  data. 

f  fo  a 

Figures  16  and  17  are  for  the  plate  below  the  boundary  layer,  at  Uoo  =  16.7  m/ s 
and  10.7  m/s.  At  both  velocities,  at  low  air  flow  rates,  the  skin  friction 
returns  to  almost  its  undisturbed  value  at  the  downstream  probe  location 
(roughly  355).  At  higher  air  flow  rates,  the  skin  friction  reductions  are 
still  significant  at  this  location. 

Figures  18  and  19  compare  the  plate  above  the  boundary  layer  with  the 
plate  below  the  boundary  layer  data  at  =  16.8  m/s  and  10.8  m/ s , 

respectively.  The  persistence  is  clearly  a  function  of  the  gravitational 
orientation  and  the  tunnel  velocity.  The  slope  of  the  curves  is  less  steep 
when  the  plate  is  above  the  boundary  layer.  Extrapolation  would  indicate 
that  the  reductions  persist  even  to  60-705.  This  increased  effectiveness 
is  in  keeping  with  other  observations. 

(e)  Upstream  Effects  due  to  Introducing  Microbubbles  in  the  Flow 

In  order  to  check  for  possible  upstream  effects  of  introducing  micro¬ 
bubbles  into  the  boundarv  layer,  an  additional  hot  film  (probe  1)  was  mounted 
10  mm  upstream  of  the  porous  injection  section.  Skin  friction  measurements  at 
this  location  are  shown  in  Fig.  20  as  a  function  of  the  air  flow  rate  for 
several  different  test  runs.  The  skin  friction  remains  unchanged  with  air 
flow,  within  the  scatter  of  the  data,  showing  that  the  introduction  of 
microbubbles  into  the  flow  is  not  felt  upstream  over  the  range  of  velocities 


presented  here. 
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( f )  Comparison  of  Hot  Film  and  Force  Balance  Measurements 

The  local  skin  friction  measurements  from  the  array  of  hot  films  can  be 
integrated  and  compared  to  the  force  balance  results  reported  in  MDM1 .  The 
integrations  were  carried  out  by  subdividing  the  measurement  section  into 
discrete  (nonuniform)  intervals  around  each  probe  location  and  assuming  the 
skin  friction  at  each  section  was  represented  by  the  measurement  at  the 
corresponding  probe.  Other  integration  procedures  are  possible  but  yield 
quantitatively  similar  results  (Madavan  1984).  The  comparisons  are  shown  in 
Fig.  21  for  different  tunnel  velocities,  plate  orientations  and  porous 
sections.  Excellent  agreement  is  noted  between  the  hot  films  and  the  force 
balance  measurements.  In  all  cases  the  maximum  discrepancy  between  the  two  is 
less  than  10%.  This  agreement  between  two  measurement  techniques  that  are 
based  on  completely  different  physical  principles  provides  substantiation  for 
both . 

(g)  Spectral  Analysis  of  Hot  Film  Signals 

The  effect  of  microbubbles  on  the  near  wall  turbulence  structure  was 
studied  by  extracting  frequency  information  from  the  hot-film  signals.  A 
typical  power  spectral  density  function,  W(f),  defined  by  the  equation 

CO 

t'  =  ;  w(f)df  ,  (5) 

w  0 

where  f  is  the  frequency  in  cycles  per  second,  is  presented  in  inner  variables 
for  no  air  injection  in  Fig.  22.  The  data  is  from  probe  2  and  the  plate  is  on 
the  bottom;  the  different  symbols  represent  various  velocities  ranging 
between  4.0  and  17.0  m/sec.  The  collapse  of  the  data  is  quite  good. 
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The  influence  of  the  microbubbles  may  be  most  strikingly  observed  by 
considering  a  sequence  of  hot-film  signal  traces  comprising  a  typical  run 
(probe  3,  16.8  m/sec,  plate  on  bottom)  as  shown  in  Fig.  23.  The  traces  have 
been  redrawn  from  the  digitized,  linearized  signals  as  processed  by  the 
computer.  The  drop  in  the  mean  level  of  skin  friction  as  the  air  flow  rate 
is  increased  from  zero  is  clearly  seen  with  reference  to  a  fixed  skin 
friction  value,  x  ^ .  In  the  presence  of  microbubbles,  the  hot-film  signals 
exhibit  a  reduced  amount  of  high  frequency  content  and  the  turbulence  energy 
appears  to  shift  toward  lower  frequencies. 

In  Figs.  24  and  25  the  power  spectral  density  functions  with  microbubbles 
for  probe  2  at  a  free  stream  speed  of  10.7  m/sec  are  shown  in  inner  variables. 
Figure  24  is  for  the  plate  below  the  boundary  layer  while  Fig.  25  is  for  the 
plate  above.  Note  that  u  and  r'  used  in  these  normalizations  are  the  values 

T  W 

with  air  injection,  so  that  the  loss  of  high  frequency  signal  is  masked  by  the 
scaling.  The  inner  variable  scaling  thus  remains  effective  for  the  micro¬ 
bubble  boundary  layer.  Similar  results  are  observed  at  higher  speeds  as  shown 
in  Fig.  26.  These  correspond  to  the  signal  traces  previously  presented  in 
Fig.  23  (probe  3,  16.8  m/sec,  plate  on  bottom). 

The  power  spectral  density  functions  indicate,  hv  the  success  of  inner 
variable  scaling,  that  the  introduction  of  air  causes  no  major  restructuring 
of  the  turbulent  boundary  layer.  We  note  that  there  is  a  strong  similarity 
between  the  results  presented  here  and  those  obtained  in  the  polymer  flow  of 
Fortuna  and  Hanratty  (1972).  Apparently  there  are  some  similarities  in  the 
drag  reduction  mechanisms  for  polymers  and  microbubbles . 
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IV.  CONCLUDING  REMARKS 

An  array  of  flush-mounted  hot  films  has  been  used  to  study  the  downstream 
evolution  and  persistence  of  the  skin  friction  reduction  in  the  microbubble¬ 
laden  turbulent  boundary  layer.  The  measurements  with  the  hot  films  are  in 
excellent  agreement  with  earlier  reported  integrated  skin  friction 
measurements  obtained  using  a  force  balance.  The  downstream  persistence  of 
the  skin  friction  reduction  beyond  the  location  of  microbubble  injection  is  a 
function  of  the  gravitational  orientation  of  the  plate  and  the  freestream 
velocity.  For  the  plate  above,  the  boundary  layer  substantial  skin  friction 
reduction  persist  for  some  355  for  low  air  flow  rates  and  up  to  60-705  for 
high  air  flows.  The  reductions  persist  for  somewhat  shorter  distances, 
roughly  506,  for  the  plate  below  the  boundary  layer. 

Contrary  to  earlier  Soviet  results,  injection  pore  size  was  found  to 
have  no  major  effect  on  the  amount  of  skin  friction  reduction.  A  hot-film 
probe  located  10  ram  upstream  of  the  porous  injection  section  indicates  that 
gas  injection  causes  no  significant  upstream  influence  for  the  range  of 
freestream  speeds  tested. 

Spectral  measurements  indicate  a  dear  loss  of  high  frequency  signal 

with  air  injection.  Inner  variable  scaling  of  the  spectra  with  u  and  t' , 

T  W 

however,  seems  as  effective  with  microbubbles  as  without.  This  indicates 
that  the  major  effect  of  the  microbubbles  in  the  turbulent  boundary  layer 
may  be  to  alter  the  local  effective  viscosity  and  density  of  the  fluid, 
thus  changing  the  local  turbulent  Reynolds  number,  rather  than  introducing 
a  major  restructuring  of  the  turbulence. 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  of  an  experimental  setup. 

Figure  2.  Schematic  snowing  hot  film  probe  locations  and  LDA  measurement 
stations  . 

Figure  3.  LDA  velocity  profile  measurements  in  the  undisturbed  boundary 
layer  presented  in  inner  variables  [Note  shifted  origins] . 

Figure  4.  LDA  velocity  profile  measurements  in  the  undisturbed  boundary 
layer  presented  in  outer  variables  [Note  shifted  origins] . 

Figure  9.  Determination  of  the  virtual  origin. 

Figure  6.  A  typical  calibration  curve  for  the  f  lu<=h-mounted  hot  film 
probes . 

Figure  7.  Hot-film  measurements  of  the  skin  friction  fluctuations  in  the 
undisturbed  boundary  layer.  The  Reynolds  number  is  based  on 
the  distance  of  the  respective  probe  from  the  virtual  origin. 

Figure  8.  Hot-film  measurements  of  the  local  skin  friction  with  micro¬ 
bubbles  as  a  function  of  the  air  flow  rate  when  the  plate  is 
below  the  boundary  layer.  U  =  16.7  m/s. 

oo 

Figure  9.  Hot-film  measurements  of  the  local  skin  friction  with  micro¬ 
bubbles  as  a  function  of  the  air  flow  rate  when  the  plate  is 
below  the  boundary  laver.  U  =  10.8  m/s. 

OO 

Figure  10.  Hot-film  measurements  of  the  local  skin  friction  with  micro¬ 
bubbles  as  a  function  of  the  air  flow  rate  when  the  plate  is 
below  the  boundary  layer.  U  =  4.6  m/ s . 


-24- 


24  August  1984 
NKM: SD: CLM: lhz 


F i guru 


Figure 


Figure 


Fi  gnri 


Figure 


F  i  g  u  r  e 


Figure 


!■'  i  guru 


11.  The  effect  of  gravity  on  the  hot-fiLm  measurements  of  the  local 
shin  friction  with  microbubbles.  Comparison  of  data  taken  with 
the  plate  above  and  below  the  boundary  layer.  U  =16.8  m/s. 

GO 

12.  The  effect  of  gravity  on  the  hot-film  measurements  of  the  local 
skin  friction  with  microbubbles.  Comparison  of  data  taken  with 
the  plate  above  and  below  the  boundary  layer.  U  =  10.7  m/s. 

CO 

IT.  The  effect  of  gravity  on  the  hot-film  measurements  of  the  local 
skin  friction  with  microbubbles.  Comparison  of  data  taken  with 
the  plate  above  and  below  the  boundary  layer.  U  =  4.6  m/s. 

oo 

Is.  The  effect  of  pore  size  variations  on  the  local  skin  friction 
measurements.  Comparison  of  data  taken  with  the  0.5  pm  and 
1 00  pm  porous  sections.  U  =  16.8  m/ s . 

OO 

15.  The  effect  of  pore  size  variations  on  the  local  skin  friction 
measurements.  Comparison  of  data  taken  with  the  0.5  pm  and 
100  pm  porous  sections.  U  =  10.7  m/s. 

OO 

16.  The  downstream  persistence  of  the  skin  friction  reduction. 

The  plate  is  below  the  boundary  Laver.  U  =  16.7  m/ s . 

30 

17.  The  downstream  persistence  of  the  skin  friction  reduction. 

The  plate  is  below  the  boundary  laver.  U  =  10.8  m/ s . 

'X> 

18.  The  downstream  persistence  of  the  skin  friction  reduction. 
Comparison  of  lata  taken  with  the  plate  above  and  below  the 
boundary  layer.  U  =  In. 8  m/s. 

19.  The  downstream  persistence  of  the  skin  friction  reduction. 

C  caparison  of  lata  taken  with  the  plate  above  and  below  the 
boundary  lav**r .  r  „  =  1  •  * . 7  m/s. 
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Figure  20. 


Figure  21. 


Figure  22. 


Figure  23. 


Figure  24. 


Figure  25. 


Hot-film  measurements  of  the  local  skin  friction  upstream  of 
microbubble  introduction  into  the  boundary  layer  [Note 
shifted  origins]  . 

Comparison  of  the  integrated  skin  friction  reductions  with 
microbubbles  measured  by  the  force  balance  with  that 
evaluated  from  the  hot  films. 

Normalized  power  spectral  density  functions  of  the  hot-film 
signals  in  the  presence  and  absence  of  microbubbles. 

A  sequence  of  linearized  hot-film  signal  traces  in  the 
presence  and  absence  of  microbubbles. 

A:  no  microbubbles; 

B:  1 0 2  air,  5.%  Ct  reduction; 

C:  14%  air,  15%  Cj  reduction; 

[):  19%  air,  28%  Ct  reduction; 


E:  26%  air,  40; 

%  Cr-  reduction;  and 

F:  36%  air,  62; 

%  Cr-  reduction. 
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Schematic  showing  hot  film  probe  locations  and  I. DA 
measu  miu'ii  t  s  t  a  t  i  ons  . 
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Figure  17.  The  downstream  persistence  of  the  skin  friction 
reduction.  The  plate  is  below  the  boundary 
layer.  if  =  10.8  m/s. 
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Figure  20.  Hot- film  measurements  of  the  local  skin 
upstream  of  m i c rohuhh 1 e  introduction  in 
boundary  layer  [Note  shifted  original. 
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Comparison  of  the  integrated  skin  friction  reductions 
vith  microbubbles  measured  by  the  force  balance  with 
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Figure 
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Normal  iced  power  spectral  density  functions  of  the 
hot-film  signals  in  the  presence  and  absence  of 
microbubbles . 
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Probe  2,  Plate  on  Bottom,  0.5pm,  U  =  16.7ra/s 


lgure  -3.  A  sequence  of  linearized  hot-film  signal  traces  in 
the  presence  and  absence  of  microbubbles . 

A.  no  microbubbles;  B:  10/1  air,  54  reduction; 
C:  14/„  air,  154  C,  reduction;  D:  194  air, 

-c.,  C f  reduction;  E;  26, 3  air,  404  C^-  reduction; 
and  F;  364  air,  624  C,-  reduction. 
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